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Abstract 

We explore the influence on the Isgur-Wise function of the colour 
electric potential between heavy and light quarks in mesons. It is 
shown that in bag models, its inclusion tends to restore light quark 
flavour symmetry relative to the MIT bag predictions, and that rela- 
tive to this model it flattens the Isgur-Wise function. Results compare 
very well with observations. 



*Work supported in part (M.S.) by the KBN grant no. 20-38-09101. We also thank the 
Norwegian Research Council for a travel grant. 
tE-mail hogaasen@vuoep6.uio. NO 
*E-mail ufsadzik@ztc386a.if.uj.edu.pl 



1 



The physics of states containing a heavy quark is the subject of much 
research. This is evident because in the hmit where the mass of the heavy 
quark grows beyond any hmit, the heavy quark can be regarded as a static 
source where the heavy quark spin effectively decouples from the interaction. 
The degrees of freedom of the light quark(s) then become the main dynamical 
variables, and the system is a fine theoretical laboratory to study properties 
of bound light quarks. It is hoped that the properties one can deduce for 
such systems is of relevance for physical systems where the heavy quark is a 
6 or a c quark. 

One of the most brilliant applications of the Heavy Quark Symmetries (see 
for example |jl|] and citation therein) is the description of semileptonic decays 
of particles containing a heavy quark in terms of only one unknown quantity 
the Isgur - Wise function (IW) [0]. The dependence of this function of light 
degrees of freedom is very complicated, and establishes an essential problem 
in all calculations that leads to numerical values of, not only branching ratios 
of semileptonic decays, but also in part of two-body decays branching ratios 
and determinations of the element \Vcb\ of Kobayashi-Maskawa matrix. 

A reliable calculation of the IW function based on first principles is beyond 
our power up to now, so a variety of different models must be used to estimate 
some numerical results. One of the possible approaches in the calculation of 
the IW function was based on the MIT bag model ||^. Bag models were 
successful in describing several important features connected with particles 
containing one heavy quark (for example in papers @], p, [0]). In this 
paper we deal with the Coulomb-like bag model, a very natural candidate 
for reasonable description of heavy - light quark system. 

In the infinite mass limit the Heavy Quark Effective Theory ([^) predicts 
that the IW function is related to the overlap of the light degree of freedom 
wavefunction multiplied by a kinematical factor. In the case of the ground 
state to ground state transitions we have 



where the ket in above equation describes the light degree of freedom wave- 
function of the parent ($q«) and the bra of the daughter ($q'J particle. The 

-'^In case of A baryons the IW function is often defined in literature also as a pure 
overlap, without kinematical factor. 




(1) 
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IW function is independent of the spin of the heavy quark. Before we come 
to our main topic we make some general remarks about the IW function vahd 
for all models satisfying the following two assumptions: 

• The valence quark approximation 

• Spherical symmetry of the ground state probability density 

Due to the first assumption the wavefunction of the light quark in the particle 
rest frame can always be written in the form: 

$0(a;) = $0(x)e-'^* (2) 

where $'^(5;) is a bispinor field and E is its ground state energy. 
The second assumption states that the product: 

p(r) = <l>°t(x)$0(a;) (3) 

is independent of angles. 

Consider the IW function for mesons. In the decay of one particle to 
another one can always choose the frame where two particles move along the 
z-axis with equal and opposite velocities (modified Breit frame 0). In this 
frame the overlap takes the form: 

< "^Q'll^Qi >= l^^d^x^^^,^ix)^Qiix)\t=o (4) 

where: 

^Qiix) = Si-v)^liix, y, jz)e-'^^^^ (5) 

is a wavefunction describing the light quark in the parent particle that is 
moving with velocity f = l^l along the negative z-axis. 

<l>Q.,{x)=S{v)<l>l,,{x,y,^z)e^^'y^^ (6) 

is a wavefunction describing the light quark of the daughter particle which 
is moving with velocity v in positive direction along the z-axis. 
Both functions are related to the wavefunction in their rest frames via boost 
operator S{v) and Lorentz transformations: 

L^liO, x) = {tivz, X, y, -fz) (7) 
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We define the volume CK as the region explored by both wavefunctions in 
the moment of the decay. In case of the MIT bag model it is the intersection 
of contracted bags in case of other models it may be the "infinite" region. 
Lets also note, that according to Heavy Quark Symmetries the wavefunctions 
$^/^ and are the same. 

If one writes down the overlap in terms of the wavefunctions from their 
rest frames, then in the integral over space there will always appear a product 
of boost matrices: {v)S{—v). Using the identity between the velocity \v\ 
and scalar product of velocities of both particles uo = vv valid in the "Breit 
frame" : 

one can write: 

S\v)S{-v) = Y.&n{^T (9) 

n=0 ^ + J- 

where only the first operator is known. It is equal to unity. Using (|l|), 
®) ©) ® s-iid ® ^iid rescaling the z - variable by a factor 7 one gets: 

e(^) = (^)[/ , d'xp{r)j,{2EJ^r) + 
oj + I Jk{o,r) V + 1 

(10) 

uj + 1 Jk{o,r) 

where -ft"(0, R) is the region explored by the light degree of freedom wave- 
function in the rest frame of particle. For the simplest bag models it is a 
spherical bag but for other models it can be an "infinite" region, jo is the 
usual spherical Bessel function of order zero. Now it is straightforward to 
find the slope parameter: 

= -^U=i = l + ]:E'<r'>-U rfWt(f)ai$°(f) (11) 
auj 2 6 2 Jk{o,r) 

Unfortunately these results are not giving us enough information. The addi- 
tional assumption that S{v) is hermitian (for example true for the free boost) 
leads to the following formula: 

d„(r) = for n=l,2,... (12) 
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and consequently: 



e(^) = (^) / , ^ d'xpir)U2EJ ^r) (13) 
00 + 1 Jk{o,r) V + 1 

p' = l + lE^<r'> (14) 

These relations were found in paper |p. The above results have a practical 
meaning for all models in which the density and the ground state energy of 
the light degrees of freedom are known. 

It was shown earlier by K. Zalewski and one of us 0, that in the case 
of mesons the MIT bag model wave functions for light quarks gave results 
for the IW function that compared very well with what we believe is known 
experimentally. When we speak of the MIT bag model here we speak of 
its simplest tractable version, the static approximation where the bag is 
spherical in its rest frame. The confinement mechanism then correspond 
to a world scalar potential which is zero inside and infinite outside the bag 
surface. There is no four vector like potential between the "infinitely" heavy 
quark, which defines the center of the bag, and the light quark. Up to 
maximal distance the bag radius, the light quark satisfies the free particle 
Dirac equation. One certainly would believe that there is an influence of the 
gluonic interaction at intermediate distances, theoretically one would expect 
something like a Coulomb-like potential, modified at very small distances by 
asymptotic freedom effects. This certainly is the result from lattice gauge 
calculations of the potential between two heavy quarks. It is therefore of 
some interest to find out how such an inner potential influences the successful 
prediction of the IW function that was obtained from the simplest (MIT) 
bag wave functions for the light quark. We therefore model the interaction 
of the heavy and light quark with a Coulomb potential for r < R, solve the 
Dirac equation there and quantize the quark energy with the MIT-Bogoliubov 



boundary condition ai r = R |]T0[. We normalize the potential so that it is 
zero at the boundary. We shall comment later on the influence of asymptotic 
freedom effects at very small r. The inner potential is now of the form: 

V{r) = -b(-- ^) (15) 



.r R 

We now define the light quark energy E as solution of the Dirac equation: 
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i/$0(r) = E$°(r) (16) 

where: 

H = dp + (3m + V{r) (17) 

Here m is the mass of the hght quark. The quark energy is determined from 
the boundary condition at the bag radius R: 

i^r^^R) = $°(i?) (18) 

where f is a radial unit vector. The hght quarks wave function and its energy 
are evidently functions of the strength of the central Coulomb potential char- 
acterized by the constant b. This is of course not a parameter that we know 
the size of. For illustrative purposes we shall use the value obtained from 
a fit to charmed mesons in 0, 6 = 0.542 and their corresponding bag radii 
for strange and nonstrange quarks. This value of b correspond to a strong 
coupling constant = 0.407, a quite reasonable value. All parameters now 
fixed, we easily compute the IW function from formula (|13|). 

In fig. 1, we show the IW function for massless light quarks in the cases 
b = (MIT bag §), b = 0.542 [|] and the ARGUS data fni. To compare the 
IW function with the data we fit the Cabbibo-Kobayashi-Maskawa matrix 



parameter \Vcb\ = 0.0386y rB/1.29ps. The of this fit is 9.3 for 7 degrees of 
freedom. As the IW functions we used the approximation of the exact result 
by the functions: 

e^(^) = i^y-''^"^ (19) 



'a; + 1' 



for B — > DeUe and 



e^i^) = (^)^-^^+^ (20) 

for Bs DsCVe similar for ones used in . The appropriate results for slope 
parameters are: 

p| = 0.98 = 1.135 (21) 

The slope parameter for B — > Deve is in a good agreement with what 
one finds in lattice calculations. One group [12] found p\ = 1.0(8), another 
13 p| = 1.2(+7-3). 
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As can be seen from the figure, the differences between the free MIT bag 
and Coulomb-like bag are not very big, but the central potential clearly has 
led to a flatter IW function. 

Considerable insight can be gained by looking at the expression for the 
derivative of the IW function at minimum recoil, i.e. at u = 1 given in 
formula (|T^), where the product X of the quark energy and the root-mean 
square radius of the light quark probability distribution enters in an essential 
manner: 

X = EV< r2 > (22) 

In the region of interest to us this product is given to a very good approxima- 
tion by a linear fit X = ai + a2mR where the coefficients ai and 02 decreases 
with increasing strength of the Coulomb potential b. For 6 = we find: 

X = 1.49 + 0.29mR (23) 

for b = 0.542: 

X = 1.15 + 0.17mR (24) 

This has interesting consequences that will show up in a moment, when we 
compute decay probabilities: The Coloumb potential tend to restore flavour 
symmetry in decay probabilities. 

Using the formulas for decay rates (for example from [Q]) and the functions 
(|0D we get: 



BriB Dlv) = 1.79(rB/1.29ps)% (25) 

BriB D*W) = 5.05(rB/1.29j9s)% (26) 

Br{Bs DJu) = 1.7A{tbJ1.29ps)% (27) 

Br{Bs DIW) = 5.02{tbJ1.29ps)% (28) 

Agreement with available data (1.7 ± 0.4%) and (4.8 ± 0.6%) [|14| for B to 



D and B to D* transitions is quite satisfactory. What is more interesting we 
also found that: 

Br{B Dlv) ^ Br{B, DsW) (29) 

Br{B -> D*liy) ^ Br{Bs DIW) (30) 

as should be expected from naive observation based on SU{3) flavour sym- 
metry. This result was absent in case of the MIT bag. It can be explained 
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by the fact that the increase of the energy of the ground state with the in- 
creasing hght quark mass is much better compensated by the decrease of the 
mean square radius in case of the Coulomb-hke than in the MIT bag. This 
is the source of the difference between (p3|) and (^41). The most important 



thing, however, is that the product X is smaller in the Coulomb bag than in 
the MIT bag so that the slope given by eq ( p!4D does not vary as much with 
varying light quark mass. 

Measuring and comparing branching ratios to test light quark flavour 
symmetry is therefore extremely interesting. 

The sceptical reader might be worried by the fact that the Coulomb po- 
tential does not incorporate asymptotic freedom, one of the most celebrated 
properties of gluonic interaction. There are all reasons to check, whether the 
singularity of the Coulomb potential at the origin distorts our results. This 
is not the case however: We integrated the Dirac equation with the potential 
V{r) of equation (|15|) for r > Tq and constant potential V^Vq) for r < Vq. 
With ro = 0.2 GeV~^ there were no differences on the 1% level for the cal- 
culated IW-function relative to the values we found with the pure Coulomb 
potential. 
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